The short-time behavior of a classical ferromagnet with double-exchange interaction 
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We investigate the critical dynamics of a classical ferromagnet on the simple cubic lattice with 
double-exchange interaction. Estimates for the dynamic critical exponents z and are obtained 
using short-time Monte Carlo simulations. We also estimate the static critical exponents v and j3 
studying the behavior of the samples at an early time. Our results are in good agreement with 
available estimates and support the assertion that this model and the classical Heisenberg model 
belong to the same universality class. 



PACS numbers: 

I. INTRODUCTION 

In the last decade many works have been devoted to 
the study of the perovskite manganites, A\- x B x MnOs, 
where A and B are rare-earth and alkaline-earth ions, 
respectively, and x is the concentration of B. The per- 
ovskite manganites are metal oxides that exhibit some 
remarkable properties. The most noticeable of them is 
colossal magnetoresistance (CMR), an extremely large 
change in the resistivity when a magnetic field is ap- 
plied in the vicinity of the critical temperature T c . Usu- 
ally these materials have been studied on the basis of 
the double-exchange theoryiSii^ (DE). In one of these 
works, Anderson e Hasegawa 2 showed that the transfer- 
ence element is proportional to cos(9/2), where 8 is the 
angle between the ionic neighbouring spins. This result 
was recently confirmed for layered manganites^. Never- 
theless, although this theory has succeeded in explain- 
ing qualitatively CMR, some authors have argued that 
it cannot alone provide a complete description of this 
phenomenon. They suggest that, in addition to the dou- 
ble exchange, a complete understanding of these materi- 
als should include strong electron correlations!, a strong 
electron-phonon interaction 8 , or coexisting phases 9 . One 
might therefore think that double-exchange alone cannot 
explain CMR in manganites^, but this remains an open 
question. What we know is that in the study of the man- 
ganites the double-exchange theory plays an important 
role, both in the study of CMR and in explaining the 
presence of a ferromagnetic state (for x « 0.3) in doped 
manganites, furnishing the basis for describing mangan- 
ites with colossal magnetoresistence. 

The mechanism of the double-exchange consists of a 
strong intrasite exchange between localized ions (core 
spin S l c ) and non- localized electrons (E g orbitals). Based 
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on this hypothesis Zener 1 succeded in explaining in 1951 
the occurrence of the metal ferromagnetic phase in man- 
ganites. Due to a strong intrasite exchange interaction 
the spin of the electron always aligns parallel to the spin 
of the ion. When an electron hops from site i to site j its 
spin must also change from being parallel to S l c to being 
parallel to S 3 C , but with an energy loss proportional to 
the cosine of half the angle between the core spins 2 . 

In Ref. [10], the critical behavior of a classical fer- 
romagnet with double-exchange interaction was studied 
and estimates for the static critical exponents is, (3 and 
7, and the critical temperature T c were derived. These 
simulations indicate that this model belongs to the uni- 
versality class of the classical Heisenberg modelii. 

In this paper, we perform short-time Monte Carlo sim- 
ulations to explore the critical dynamics of a classical 
ferromagnet with double-exchange interaction. We eval- 
uate the dynamic exponents z, and the static expo- 
nents v and p. Our estimates for static indices are in 
good agreement with previous results available in the 
literatureifi^2iii for the same model. In addition, the re- 
sult for the dynamic exponent z compare well with pre- 
vious estimates for the classical Heisenberg model. In 
summary, present results confirm that both models be- 
long to the same universality class. 

The article is organized as follows. In the next Section 
we define the model. In Section III we give a brief descrip- 
tion of short-time Monte Carlo technique and present our 
results; in Section IV we summarize and conclude. 



II. THE MODEL 

We use the short-time Monte Carlo simulations to in- 
vestigate the dynamical critical behavior of a classical 
spin model with double-exchange interaction. We con- 
sider a simple cubic lattice Lx Lx L with periodic bound- 
ary conditions. The Hamiltonian for such system is given 
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n = ~Y J A/ 1 + s * -S,-, (i) 

<i,3> 

where < i,j > indicates that the sum runs over all 
nearest-neighbor pairs of lattice site, J is the ferromag- 
netic coupling constant and the spin Si = (Sf , Sj, Sj£) is 
a three-dimensional vector of unit length. 

We used lattice sizes L = 20, 25, 30, 40 and 60 
and performed short-time simulations at the critical 
temperature^ 

T c = 0.74515, 

in units of J/ks, where ks is the Boltzman's constant. 

Our estimates for each exponent were obtained from 
five independent bins. The dynamic evolution of the 
spins is local and updated by the Metropolis algorithm. 
For 20 > L > 40 the estimate for the exponent 8 was per- 
formed using 10000 independent samples, while for the 
exponents z, v and (3 we used 5000 samples. For L = 60 
all the estimates were performed using 2000 samples. In 
order to maximize the quality of the linear fitting, each 
exponent is determined in the time interval [200,500], 
except for the exponent z which was measured in the 
interval [30,250]. 

III. SHORT-TIME SCALING RELATIONS AND 
RESULTS 

Following the works of Janssen, Schaub and 
Schmittmanni^ and Husei^, the study of the critical 
properties of statistical systems became in some sense 
simpler, because they allow to circumvent the well-known 
problem of the critical slowing down, characteristic of 
the long-time regime. They discovered using renormal- 
ization group techniques and numerical calculations, re- 
spectively, that there is universality and scaling behavior 
even at the early stage of the time evolution of dynami- 
cal systems without conserved quantities (model A in the 
terminology of Halperin et al—) . 

In this new universal regime, in addition to the fa- 
miliar set of static critical exponents and the dynamic 
critical exponent z, a new dynamic critical exponent 8 is 
found. This new critical index, independent of the previ- 
ously known exponents, characterizes the so-called critical 
initial slip, the anomalous increase of the magnetization 
when the system is quenched to the critical temperature 

Several works on phase transitions and critical phe- 
nomena using this technique have been published^. The 
results corroborate the prediction of universality and 
scaling behavior in the short-time regime and the static 
critical exponents so obtained are in good agreement with 
well-known values calculated in the equilibrium. In addi- 
tion, this approach has permitted obtaining more precise 
estimates for the dynamic critical exponent z. 



The dynamic scaling relation obtained by Janssen et 
al for the fc-th moment of the magnetization, extended 
to systems of finite sizei^, is written as 

M W (t, r, L, m ) = b~ kf3/v M {k ^ (b- z t, b 1/u r, b~ l L, b x °m ), 

(2) 

where t is time, b is an arbitrary spatial rescaling fac- 
tor, r = (T — T c ) /T c is the reduced temperature and L 
is the linear size of the lattice. The exponents (3 and 
v are the equilibrium critical exponents associated with 
the order parameter and the correlation length, and z is 
the dynamic exponent characterizing time correlations in 
equilibrium. For a large lattice size L and small initial 
magnetization mo at the critical temperature (t — 0), the 
magnetization is governed by a new dynamic exponent 8, 

M{t) ~ m t e (3) 

if we choose the scaling factor b = t x l z . 

In addition, a new critical exponent xq, which repre- 
sents the anomalous dimension of the initial magnetiza- 
tion mo, is introduced to describe the dependence of the 
scaling behavior on the initial conditions. This exponent 
is related to 9 as xq = 8z + f3/v. In the following sec- 
tion we evaluate the critical exponents 8 7 z, v e (3 using 
scaling relations obtained from Eq. (J2J. 

A. The dynamic exponent 6 

Usually the dynamic critical exponent 9 is calculated 
using Eq. J3J or through the autocorrelation 

A(t)~<M, (4) 

where d is the dimension of the system. 

In the first case, careful preparation of the initial con- 
figurations (mo <C 1) is essential, as well as the delicate 
limit mo — > 0. In the autocorrelation case, one must to 
know in advance the exponent z, which is an order of 
magnitude greater than 8, so that a small relative error 
in z induces a large error in 8. 

In the present work we find the dynamic critical expo- 
nent 9 using the time correlation of the magnetization 19 

C(t) = (M(O)M(i)) ~ t e . (5) 

Here, the averaging is over a set of random initial con- 
figurations. The time correlation allows the direct calcu- 
lation of the dynamic exponent 8, without the need of a 
careful preparation of the initial state nor of the limiting 
procedure, the only requirement being that {Mo} = 0. 

In Fig. ^ we show the time dependence of the time 
correlation C(t) in double-log scale for the system with 
L = 60. Table IJ shows our estimates for the critical 
exponent 8 along with d/z, 1/vz and (3/vz. 

These results show that for L > 25 the finite size ef- 
fects are less than the statistical errors. Therefore, we 
conclude that the values for an infinite lattice are within 
the errorbars of our results for L = 60. 
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FIG. 1: Time evolution of the time correlation of the total 
magnetization on log scales. Errorbars are smaller than the 
symbols. Each point represents an average over 5 sets of 2000 
samples. 



FIG. 2: The cumulant U 4 (t, L) for L = 40 and 60 plotted 
versus time t with initial states mo = 1. The line shows the 
lattice L — 60 re-scaled in time. The cumulant for both lattice 
sizes was obtained from five independent bins. For L = 40 we 
used 5000 samples and for L — 60 we used 2000 samples. 



TABLE I: Critical exponents for the double-exchange model. 



d/z 1/vz P/vz 



20 0.460(4) 1.499(1) 0.766(6) 0.276(3) 

25 0.473(4) 1.515(2) 0.734(5) 0.267(3) 

30 0.479(8) 1.521(5) 0.748(7) 0.267(2) 

40 0.482(4) 1.525(7) 0.743(5) 0.264(1) 

60 0.480(7) 1.519(8) 0.738(7) 0.263(2) 

B. The dynamic exponent z 

In the short-time regime the dynamic critical exponent 
z can be evaluated using various methods, such as the 
autocorrelation 20 (Eq. (0}, with the previous knowledge 
of 9), or the second moment M^ 2 \t)r^ (when the static 
exponents j3 and v are known). In the present work we 
estimate the exponent z in two independent ways, us- 
ing two different techniques. First, we used the time- 
dependent fourth-order Binder cumulantiSiSi 



Ui(t,L) 



1 



(M 4 ) 
3(M 2 } 2 



(6) 



with ordered initial configurations (mo = 1). If r = 
this equation depends only on t/L z , according to scaling 
laws valid in the beginning of the evolution. 

In this case, according to the scaling relation Eq. 
the cumulants obtained in two different lattices obey the 
equation 



U 4 (t,L 1 ) = U i (b-%L 2 ) 



(7) 



with b = L\j L 2 . This equation shows that the exponent 
z is obtained by finding the time re-scaling factor b~ z 
which colapses curves for two lattices of different sizes. 

In Fig. |2jwe show the Binder cumulant as a function of 
the time for the lattices L\ = 60 and L 2 = 40. The line 



shows the colapse of the larger lattice re-scaled in time. 
Our estimate for z, obtained through the x 2 test2&, for 
Li = 60 and L 2 = 40 is 



z = 1.98(3) 



(8) 



This result is consistent with that obtained by Peczak 
and Landau^ for the classical Heisenberg model. Using 
equilibrium tecniques, these authors found z — 1.96(7). 
We also estimate the exponent z for other lattice sizes 
and our results are consistent with the Eq. JSJ. For 
L\ = 40 and L 2 = 25 we obtained z = 1.96(2), whereas 
for L\ = 60 and L 2 = 30 we obtained z = 1.96(4). 

A more precise estimate for the dynamic exponent 
z can be obtained by combining results from samples 
submitted to different initial conditions (mo = 1 and 
mo = 0). According to the Ref. [24] the function 
F 2 (t) = M^(t)/(M(t)) 2 behaves as t d / z where d is 
the dimension of the system. This approach proved to 
be very efficient in estimating the exponent z, accord- 
ing to results for the Ising model, the three- and four- 
state Potts models^, the tricritical point of the Blumc- 
Capel model 2 ^, the Baxter- Wu models and nonequilib- 
rium models like Domany-Kinzel and contact process 2 ^. 
In this technique, for different lattice sizes, the double- 
log curves of F 2 versus t fall on a single straight line, 
without any rescaling of time, resulting in more precise 
estimates for z. 

The time evolution of F 2 is shown on log scales in Fig. 
[3]for L = 60. In Table^we show the values of d/z for the 
five lattice sizes. Taking into account the value of this 
ratio for the lattice L = 60 we obtain 



1.975(4) 



(9) 



Our results for z indicate that besides presenting the 
same static critical exponents double-exchange and 
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FIG. 3: Time Evolution of Fi- The error bars, calculated over 
5 sets of 2000 samples, are smaller than the symbols. 



Heisenberg models have the same dynamic exponent 
when submitted to the same updating. 



C. The static exponents v and 

The static exponent v can be obtained by fixing b~ z t 
1 in Eq. J3J and differentiating InM (t, r) with respect to 
r in the critical point. Taking into account samples with 
ordered initial configurations (mo = 1), we obtain the 
following power law 



d T lnM{t,T)\ T=0 ~t x l vz . 



(10) 



In numerical simulations we approximate the derivative 
by a finite difference; our results were obtained using 
finite differences of T c ± 5 with i5 = 0.01, in units of 
J/ks- In Fig. 0| the power law increase of Eq. I|10|) is 
plotted in double-log scale for L = 60. 

From the slope of the curve one can estimate the crit- 
ical exponent \jvz (see Table HJ. 

Using the exponent z obtained from U^t^L) and the 
estimate of 1/vz for L = 60, we obtain the static expo- 
nent 



v = 0.68(2). 



(11) 



while following the scaling relation F% (t) the result is 

^ = 0.682(8). (12) 

Finally we can evaluate the static exponent /3 by the 
dynamic scaling law for the magnetization 



M(t) ~t-P/"*, 



(13) 



obtained from Eq. (|2J), by considering large systems and 
setting b — t x / z at the critical temperature r = 0. In 
Fig. we show the time evolution of the magnetization 
in double- log scale for L — 60. The estimates of the 
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FIG. 4: The time evolution of the derivative d T \nM(t, r)\ T= o 
on log scales in a dynamic process starting from an ordered 
state (mo = 1). Error bars are smaller than the symbols. 
Each point represents an average over 5 sets of 5000 samples. 



J_ 



0.263(2) 
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FIG. 5: The time evolution of the magnetization for initially 
ordered samples (mo = 1). The errorbars calculated over 5 
sets of 2000 samples, are smaller than the symbols. 



exponent (3/vz for the five lattice sizes are shown in Table 
III 

Using the previous result obtained for 1/vz (Table QJ, 
we find 



= 0.354(5). 



(14) 



Our estimates for the static exponents v and j3 are in 
good agreement with the theoretical results of the Ref. 
[10] {v = 0.6949(38) and = 0.3535(30)), Ref. [11] 
(u = 0.7036(23) and (3 = 0.3616(31)) and Ref. [28] (y = 
0.704(6) and (3 = 0.362(4)), as well as the experimental 
results of the Ref. [29] (/? = 0.37(4)) and the Ref. [30] 
{P = 0.374(6)). 
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IV. DISCUSSION AND CONCLUSIONS 

In this article we have performed short-time Monte 
Carlo simulations in order to evaluate the dynamic and 
static critical exponents of a classical ferromagnet with 
double-exchange interaction. The exponent 9 was esti- 
mated using the time correlation of the magnetization. 
We found the dynamic exponent z through the colapse 
of the fourth-order time-dependent Binder cumulant and 
also, alternatively, using the function F 2 (t) which com- 
bines simulations performed with different initial condi- 
tions. Using scaling relations for the magnetization and 
its derivative with respect to the temperature at T c , we 



have obtained the static exponents v and (5 which are 
in good agreement with the values available in litera- 
ture. Comparison of our results with those for the classi- 
cal Hcisenberg model corroborate the assertion that the 
models belong to the same universality class. 
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